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ABSTRACT 
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u 

Q I The accuracy and robustness of a simple method to estimate the total mass profile 

^ . of a galaxy is tested using a sample of 65 cosmological zoom-simulations of individual 

C/3 ' galaxies. The method only requires information on the optical surface brightness and 

^ ^ I the projected velocity dispersion profiles and therefore can be applied even in case 

of poor observational data. In the simulated sample massive galaxies (cr ~ 200 — 400 
km s~ ) at redshift z = have almost isothermal rotation curves for broad range 
J> . of radii (RMS ~ 5% for the circular speed deviations from a constant value over 

' 0.5i?cff < r < 3Rcs). For such galaxies the method recovers the unbiased value of 

^^O I the circular speed. The sample averaged deviation from the true circular speed is less 

. than ~ 1% with the scatter of ~ 5 — 8% (RMS) up to i? ~ 5i?off- Circular speed 

estimates of massive non-rotating simulated galaxies at higher redshifts (z = 1 and 
■rj" , z = 2) are also almost unbiased and with the same scatter. For the least massive 

' galaxies in the sample (cr < 150 km s~^) at z = the RMS deviation is ~ 7 — 9% and 

the mean deviation is biased low by about 1 — 2%. We also derive the circular velocity 
profile from the hydrostatic equilibrium (HE) equation for hot gas in the simulated 
galaxies. The accuracy of this estimate is about RMS ~ 4 — 5% for massive objects 
(M > 6.5 X and the HE estimate is biased low by ~ 3 - 4%, which can be 

traced to the presence of gas motions. This implies that the simple mass estimate can 
^ I be used to determine the mass of observed massive elliptical galaxies to an accuracy 

■ - - of 5 — 8% and can be very useful for galaxy surveys. 

Key words: Galaxies: Kinematics and Dynamics, X-Rays: Galaxies 



1 INTRODUCTION 

The accurate determination of galaxy masses is a crucial 
issue for galaxy formation and evolution models. Disentan- 
gling dark matter and baryonic matter of a galaxy permits 
testing the predictions of ACDM-cosmology and probing the 
mass function. An algorithm for deriving the mass of a spi- 
ral galaxy is straight forward - one just need to measure 
a rotation curve from gas or stars that can be safely as- 
sumed to be on circular orbits. For elliptical galaxies the 
situation is less simple. There is no 'perfect' (in terms of ac- 



curacy) tracer to measure the total gravitational potential. 
The main problem is the degeneracy between the anisotropy 
of stellar orbits and the mass. The shape of stellar orbits 
is not known a priory and different combinations of or- 
bits may give the same distribution of light. Several differ- 
ent approaches for mass determination were proposed and 
succesfuUy implemented, like strong and wea k lensing (e.g. 
iGavazzi et al.1120071 : iMandelbaum et al.ll2006l'). modelling of 



X-ray emission of hot gas in galaxies (e.g. Humphrey et al.l 
l200d : IChurazov et al.ll2008l '). Schwarzschild modelling of stef 
lar orbits, etc. Accurate data on the projected line-of-sight 
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velocity distribution with information on higher-order mo- 
ments enables an accurate determina tion of the ma s s dis- 
tribution for nearb y ellipticals (e.g. iGerhard et al. 1 1 19981 : 
I Thomas et aLllioilh . However, in case of minimal available 
data detailed modelling is often not possible. Therefore it is 
important to find a method to measure galaxy masses with 
reasonable accuracy which gives an unbiased estimate when 
averaged over a large number of galaxies. In particular, it 
can be extremely useful while analysing large surveys, espe- 
cially at high redshifts when detailed observational data of 
each individual galaxies are often not available. 

The simplest way of estimating the mass of a galaxy is 
base d on the projected velo city dispersion in a fixed aperture 
fe.g. lCappellari et al .120061). A slightlv more co mplicated ap- 
proach is described in Ichurazov et al.l (|2010l ). To estimate 
the mass the only information required is the light profile 
and either the dispersion profile measurement or at least a 
reliable dispersion measurement at some radius. Testing this 
particular method on a sample of simulated galaxies is the 
subject of this paper. The main questions that we want to 
address are (i) What is the accuracy of this method? (ii) 
Does it give an unbiased result? (iii) What are the restric- 
tions for application of this method? 

The structure of the paper is as follows. In section[21 we 
provide a brief description of the method. In section |3] we 
describe the sample of simulated galaxies which is used to 
test the method. The analysis of the accuracy of the method 
is presented in section [4] where we also discuss alternative 
methods for determining the circular velocity. A summary 
on the bias and accuracy of the various methods is given in 
section [5] with conclusions in section [S] 



2 DESCRIPTION OF THE METHOD 



The m ain idea of the method is described in lChurazov et al] 
Here we just provide a brief summary. 
The method is based on the stationary non-streaming 
spherical Jeans equation: 



or r dr 



(1) 



where j(r|3 is the stellar luminosity density, cjr{r) is the 
radial component of the velocity dispersion tensor (weighted 
by luminosity), /3(r) =1 — og/a^ \& the stellar anisotropy 
parameter [ae — cr^ because of the assumed spherical sym- 
metry) and $(r) is the gravitational potential of a galaxy. 

While the stellar luminosity density j(r) and radial dis- 
persion ar{r) can not be observed directly they contribute 
to the two-dimensional surface brightness I{R) and the ve- 
locity dispersion <j{R) profiles: 
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^ Throughout this paper we denote a projected 2D radius as R 
and a 3D radius as r. 



Assuming /3(r) = const we note that /3 = for systems 
where the distribution of stellar orbits is isotropic, /3 = 1 if 
all stellar orbits are radial and /3 — >■ — oo if the orbits are 
circular. 

Assuming the logarithmic form of the gravitational po- 
tential $(r) = ln(r) + const and using local properties 
of given /(i?) and cr(ii) one can calculate a circular veloc- 
ity Vc for three different types of stellar orbits: isotropic 
(or — — ae, 13 = 0), radial [a^ — ae = 0, /S = 1) and 
circular (cr^ = 0, /3 — > — oo). These relations are given by: 

Vr = ai,oiR) ■ v/l + Q + 7 



K""='^circ(i?)M/2 



1 + a + 'y 



(4) 



where 

dlna 



^ d In J(i^) 
"~ dlnR ' 



dlnR 



2 ^_ d^\n[I{R)a^] ^ 



d(lni?)2 



In case of noisy data on the dispersion velocity profile 
the subdominant terms 7 and S can be neglected, i.e. the 
dispersion profile is assumed to be fiat, and equations Q 
are simplified to: 



14"° = aieo(fl) ■ V^TT 



a -I- 1 



(6) 



Let us call a sweet spot the radius at which all three 
curves V^^°{R), V;:"''{R) and V^'^'^iR) are very close to each 
other. One can hope that at the sweet spot the sensitiv- 
ity of the method to the stellar anisotropy parameter /3 is 
minimal and the estimation of the circular speed at this par- 
ticular point is reasonable. E.g. from equations ((6]) it is clear 
that in case of the power-law surface brightness profile with 
a = 2 and /? — const the relation between the circular speed 
and the projected velocity disp ersion does no t depend on 
the anisotropy parameter (e.g. iGerhardI 1 19931 ). While the 
derivation of equations Q, ((6]) relies on the assumption 
about a fiat circular velocity profile, tests on model galaxies 
with non-logarithmic potentials, non-power law behaviour 
of the surface brightness and line-of-sight velocity disper- 
sion profiles and with the anisotr opy parameter P varying 
with radius (|Churazov et al.ll2010l ) have shown that the cir- 
cular speed can still be recovered to a reasonable accuracy. 
Now we extend these tests to a sample of simulated elliptical 
galaxies. 

This method for evaluating the circular speed is not 
only simple and fast in implementation but it also does 
not require any assumptions on the radial distribution of 
anisotropy /3(r) and mass M{r). 

The mathematical de r ivatio n of equations (|4I6|I can be 
found in IChurazov et al.l (|2010l ). A similar approach and 
analytic formulae for kinematic de projection a n d ma ss in- 
version also can also be found in IWolf et all (|2010| j and 
iMamon et al] (|2010l '). 
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3 THE SAMPLE OF SIMULATED GALAXIES 
3.1 Description of the sample 

Simulations provide a useful opportunity to test different 
methods and procedures as all intrinsic properties of a sys- 
tem at hand are known. The main drawback of simulated 
objects is that they may not include all physical processes 
that take place in reality and thus may not reflect all com- 
plexity of nature. To test the procedure under considera- 
tion we have used a sample of 65 cosmologic al zoom sim- 
ulations partly presented in lOser et all (|201(]| '). These SPH 
simulations include feedback from supernovae type II, a uni- 
form UV- background radiation field, star formation and ra- 
diative Hydrogen and Helium cooling but do not include 
ejective feedback in the form of supernovae driven winds. 
Present-day stellar masses of simulated galaxies range from 
2.18 X lOi°M0/i-i to 28.68 x lQ^°MQh-^ inside 30 kpc. 
The softening length used in simulations is about Rsoft=4:00 
pc , h = 0.72. Typically the softening can affect profiles 
up to ~ SRsoft, which is ~ 1.7 kpc in our case. We have 
followed a conservative approach and restricted the analysis 
to radii larger than 3 kpc. It should be noted that low-mass 
simulated galaxies may have no real counterparts possibly 
due to lack of important physical processes (e.g., significant 
winds) in s i mulat ions. However, it has been demonstrated in 
lOser et ahl (|201ll ) that the massive simulated galaxies have 
properties very similar to observed early-type galaxies (see 
also Figure 4), i.e. they follow the observed scaling relations 
and their evolution with redshift. For detailed descripti on of 
simulations and included physics see lOser et al.l |20l3)- 

To effectively increase the number of galaxies we have 
considered three independent projections of each galaxy. So 
the whole sample of simulated galaxies consists of 195 ob- 
jectfl 



3.2 Isothermality of potentials in massive galaxies 

First of all we have found that massive galaxies in the sam- 
ple have almost isothermal rotation curves over broad range 
of radii. To demostrate this statement (Figure [ij we have 
selected galaxies with a projected velocity dispersion at the 
effective radius CT(-Rofi) (procedure of computation _Rcff is 
described in section [3^ greater than 200 km s~^ and plot- 
ted their circular velocity curves Vc = \/ GM{< r) /r as a 
function of r/Rcs- G is the gravitational constant, M{< r) 
is the mass enclosed within r and RcS is the effective radius 
of the galaxy. The circular velocity curves were normalised 
to the value of Vc averaged over r £ [0.57ieff, 37ieff]- Three 
circular velocity curves that make the most significant con- 
tribution to the RMS actually correspond to galaxies with 
the effective radius R^ff < 6 kpc. The fact that for these 
galaxies 0.5iieff is close to the softening length may affect 
the scatter. 



^ Nevertlieless, for calculating an error in a bias estimation (= 
RMS /^/N) we conservately use the number of galaxies rather 
than the number of projections as the subsamples corresponding 
to dilTerent projections are not entirely independent. 




Figure 1. Circular velocity curves of massive galaxies (cr(_Rj,fj) > 
200 km s^^) as a function of radius r. Individual rotation curves 
normalised to the speed averaged over [0.5-Roffi 3-Rcft] are shown 
in black, green dashed lines indicate the interval [1 — RM 5, 1 -|- 
RMS], where RMS = 4.9%, the red thick line represents the 
overall trend Vc oc r^^-Oe 



3.3 Analysis procedure 

The analysis of each galaxy consists of several steps de- 
scribed below. 

Step 1: Excluding satellites from the galaxy image. 

Usually an image of a simulated galaxy (the distribu- 
tion of stars projected onto a plane) contains many satel- 
lite objects and needs to be cleaned. Exclusion of satellites 
makes the surface brightness and the line-of-sight velocity 
dispersion profiles smoother and reduces the Poisson noise 
associated with satellites. The algorithm we used for remov- 
ing satellites is as follows: first, for each star a quantity w 
characterising the local density of stars {w oc p« ^''^) and 
analogous to the HSML (the SPH smoothing length) was 
calculated and the array of these values was sorted. Then the 
(0.4- A^stars)"' term of the sorted lo-array was chosen as a ref- 
erence value Wo- Natara is a total number of stars in a galaxy 
and a factor in front of Natara is some arbitrary parameter 
(the value 0.4 was chosen by a trial-and-error method). Stars 
with the 3D-radius r > 10 kpc and w < Wo are considered as 
members of a satellite. After projecting stars onto the plane 
perpendicular to the line of sight we have excluded all satel- 
lites together with an adjacent area of 1.5 kpc in size. The 
inititial and final images of some arbitrarely chosen galaxy 
(the virial halo mass is ~ 1.7 x W^"^ MQh~^) are shown in 
Figure [3 

Step 2: Evaluating I{R) and a{R). 

All radial profiles have been computed in a set of loga- 
rithmic concentric annuli around the halo center. To calcu- 
late the surface brightness profile, corrected for the contam- 
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Figure 2. Excluding the satellites. 150 kpc X 150 kpc. Left: Initial galaxy image. Right: Cleaned galaxy image. 




Figure 3. Influence of satellites on the surface brightness (the upper panel) and the projected velocity dispersion profiles (in the middle). 
Open black circles correspond to the initial galaxy image and solid red circles - to the galaxy image without satellites. The black dashed 
curve is the smoothed curve for the initial data and the black solid curve is for the cleaned data. The bottom panel shows the true 
circular velocity (black thick line) and recovered circular velocity for the isotropic distribution of stellar orbits (in blue) for initial data 
(dashed) and cleaned data (solid). It is clear that removing satellites reduces the scatter in the line-of-sight velocity dispersion data and 
makes the profile smoother. 



ination from the satellites, we have first counted the number 
of stars in each annulus, excising the regions around satel- 
lites. The surface area of each annuli has been also calcu- 
lated, excluding the same regions. The ratio of there quanti- 
ties gives us the desired 'cleaned' surface brightness profile. 
The average line-of-sight velocity of stars and the projected 
velocity dispersion have been calculated similarly. 

Importance of the 'cleaning' procedure and the resulting 
profiles of I{R) and (j{R) are shown in Figure[3] The surface 
brightness data (open circles correspond to the initial ('un- 
cleaned') image and red solid circles to the 'cleaned' image) 



and the smoothed curves (the calculation of these curves is 
described in Step 3) are shown in the upper panel, the pro- 
jected velocity dispersion profiles are shown in the middle 
panel. The true circular velocity Vj'^"'^(r-) (black solid curve) 
and recovered from the initial data (blue dashed line) and 
from 'cleaned' data (blue solid line) circular velocity for the 
isotropic distribution of stellar orbits V^^° (the first equation 
in (m) are shown in the bottom panel. The last curve is in 
better agreement with the true velocity profile. All results 
and figures in this paper are restricted to the region R > 3.0 
kpc. 
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Step 3: Taking derivatives. 

To take deriv a tives we follow the procedure described 
in lChurazov et al] (|2010l ) in Appendix B. The main idea is 
that all data points participate in calculating the derivative 
but with different weights. The weight function is given by 



1.5 



W{Ro,R) = exp 



{InRo-lnRY 



2A2 



(7) 



where 7?o is the radius at which the derivative is being cal- 
culated and the parameter A is the width of the weight 
function. 

Both observed and simulated surface brightness profiles 
are typically quite smooth so we have used A/ — 0.3 to 
calculate the logarithmic derivative ciln/(_R)/dln_R. For the 
line-of-sight velocity dispersion data we have used A^ — 0.5. 
With the assumed values of A the local perturbations are 
smoothed out but the global trend of the profiles is not af- 
fected. Changing values A/ and Ao- in the range [0.3, 0.5] 
does not significantly influence our final resul10. The differ- 
ence (in terms of circular velocity) is less than 1%. As an 
example the smooothed curves for the I{R) and (j{R) data 
in Figure [3] are calculated using this procedure. 

We have also tested the influence of parameters of the 
presented smoothing algorithm. As long as the smoothed 
curve describes data reasonably well neither the functional 
form of the weight function nor other parameters (like higher 
order terms in expansion Inl(R) = a(lnf?)^ -1- b\nR -|- c or 
cr(_R) = a{lnR)^ + bin R + c) significantly affect the final 
result. 

Step 4: Estimating the circular velocity. 

Applying equations Q or ([6]) to the smoothed I{R) and 
cr(i?) we have calculated Vc-profiles assuming isotropic, ra- 
dial and circular orbits of stars. Then we have found a radius 
(a sweet point i?sweet) at which the quantity (14'°° — V)^ + 
{Vr'i-Vf + {Vr''-Vf, whereF = (V,'=°-FKf"''+K"")/3, 
is minimal. The value of the isotropic velocity profile at this 
particular point is the estimation of the circular velocity 
speed we are looking for. We take V^^° as an estimate of the 
Vc{R) (rather than Vf'"' or V^'"^) for two reasons. Firstly, at 
around one effective radius the dominant an isotropy for most 
ellipti cal galaxies is a^z < ctrr ~ cr^^ (jCappellari et al.l 
(|2007l )). The spherically averaged anis o tropy is therefore 
only moderate (see also ICerhard et al.l (|200fl ). Figure 4). 
Massive elliptical galaxies are the most isotropic. Thus an 
isotropic orbit distribution is a much better approximation 
than purely radial or circular orbits. Secondly, the value of 
Vc"^° is less prone to spurious wiggles in I{R) and a{R). 

The effective radius R^b is calculated as a radius of the 
circle which contains half of the projected stellar mass, tak- 
ing into account effects of cleaning. We found that in the 
simulated data-set the value of the effective radius depends 
on the maximal radius used to calculate the total number 
of stars in a galaxy. The problem is especially severe for 
the most massive galaxies as they have an almost power- 
law 3D stellar density distribution p, oc r~° with a ~ 3. 



^ If, however, we choose a width of the weight function smaller 
that A = 0.3 the local scatter in the data is not smoothed out 
and the results become ambiguous. 



0.5 - 



Auger et al'lO 
this work 
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Figure 4. R^g — relation. The blue solid line is the linear fit 
to data points from the simulatio ns. The green dash ed line is the 
observed mass-size relation from l| Auger et al.ll201(t . 



In our analysis (in contrast to lOser et al.l (|201lh ) we have 
not introduced any artificial cut-off and used all stars in the 
smooth stellar component (excluding substructure) of the 
main galaxies out to their virial radii for the calculation of 
the effective radius. The resulting effective radii as a func- 
tion of total stellar mass (in logarithmic scale) are shown in 
Figure Q. The slope and the normalization of t he Res — M<. 
relatio n are close to the fit of SLACS data bv lAuger et al.l 
(|20IOl ). 

The axis ratio q of each projection of a galax;y is cal- 
culated as a square root of eigenvalues of the diagonalised 
inertia tensor. The inertia tensor is computed within the 
effective radius without excluding substructures. We have 
found that q is not sensitive to our cleaning procedure as 
normally there are almost no satellites within -Rcff- 



4 ANALYSIS OF THE SAMPLE 
4.1 At a svi^eet point 

For each galaxy in the sample we have performed all steps 
described above and we have selected the radius at which 
the circular velocity curves for isotropic, circular and radial 
orbits (equations Q) intersect or lie close to each other. 
Then we have calculated the value of the isotropic speed Vc°° 
at this radius. To measure the accuracy of our estimates let 
us introduce a deviation from the true circular speed Aopt = 
(y^^° - V?-^"") /V^''"', where and V?"^"" should be taken 
at the sweet spot -Rswcct- The subscript 'opt' (= optical) 
is used to distinguish this method (based on optical data) 
from circular speed calculations based on X-ray data. We 
have plotted the number of galaxies (normalised to the total 
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Full analysis, subsample 'A' 




-20 
deviation A, 



-(V — V'-™«)/V"''", % 



Figure 5. The fraction of galaxies (in %) as a function of deviation A^pt = 

(V;iso _ t/truc-) /v;"^"<= evaluated via equations Q at different 
radii: iJswoct (panel (A)), i?eff (panel (B)), 0.5/?efi (panel (C)) and 2R^ff (panel (D)). 
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Figure 6. Left: Example of the galaxy that perfectly suits for the analysis. The surface brightness and the projected velocity dispersion 
profiles are shown in panels (A) and (B) correspondingly. Data are represented as red points and smoothed curves that were used to 
compute derivatives (a, 7, 5) as black solid lines. The auxilary coefficients o, 7, —5 and a -I- 7 are shown in panel (C) in red solid, blue 
dotted, green dash-dotted and black dashed lines, respectively. Circular velocity profiles for isotropic orbits of stars (blue solid line), pure 
radial (green dash-dotted) and pure circular (magenta dashed) orbits as well as the true circular speed (black thick curve) are presented 
in panel (D) for the full version of the analysis (equations And the same curves for the simplified analysis (equations ^) are 

shown in panel (E). Right: Example of the galaxy with large deviation Aopt due to merger activity. The crest in the projected velocity 
dispersion profile at ij ~ 20 kpc loads to the significatly overestimated value of the circular speed. 
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Figure 7. Left: Shown in cyan is the histogram for deviations for galaxies with the axis ratio q < 0.6, in black is the histogram for 
the same galaxies but seen in a projection with the axis ratio q close to unity (= seen along the rotation axis). Right: The histogram 
for deviations for the sample when merging and oblate galaxies seen along the rotation axis are excluded (subsample 'G'). The average 
deviation A^^ = (-1.2 ± 0.9)%, RMS = 6.8%. 



number of galaxies and expressed in %) versus the deviation 
Aopt in a form of a histogram. To have an idea whether 
the method under consideration gives resonable accuracy, 
histograms for deviations at i?efi, O.STieff and 27?efF are also 
shown. The whole sample ('subsample A') is presented in 
Figure [5] The sample averaged value of the deviation Aopt 
is slightly less than zero in all cases. For example, at the 
sweet point A^^ = (-1.8 ± 1.1)% while the RMS = 8.6?^. 

Large deviations (~ 30 — 40%) are seen only in galaxies 
with ongoing merger activity. The influence of mergers ap- 
pears as 'waves' in the projected velocity dispersion profile. 
The example of such a system is shown in Figure |6] (right 
panel). The presence of such 'waves' indicates that the cir- 
cular speed could be significantly overestimated (by a factor 
of ~ 1.2 — 1.5), which is not surprising as the method is 
based on the spherical Jeans equations and the assumption 
about dynamical equilibrium is violated. When the profiles 
I{R) and (j{R) are smooth and monotonic the circular speed 
can be recovered with much higher accuracy (Figure |6] left 
panel) . 

The sample includes galaxies with different values of el- 
lipticity. The axis ratio q (computed from the diagonalized 
inertia tensor within Res) ranges from 0.19 to 0.99. To test 
the possible influence of the ellipticity on the accuracy of 
estimates we have selected galaxies with axis ratio q < 0.6. 
The resulting distribution as a function of the circular speed 
deviations is almost symmetric, unbiased, with RMS ~ 8% 



(Figure[71). On the other hand, if we consider the same galax- 
ies seen in a projection with the maximum value of the axis 
ratio q, we get the distribution appreciably biased toward 
negative values of the deviation (Aopt = (—10.2 ± 1.6)%). 
The reason for this bias is rotation. When observing a galaxy 
along its rotation axis the projected velocity dispersion is 
appreciably smaller than for perpendicular directions. To 
further test this statement we have rotated each galaxy so 
that the principal axes of the galaxy {A ^ B ^ C) coincide 
with the coordinate system {x, y and z, correspondingly) 
and analysed velocity maps for each projection. As a cri- 
teria for rotation we have used the anisotropy-parameter 



{v/ay 



where v is the average rotation ve- 



N ' 



RMS 



E (^ - 
N - 1 



\/(l -J)/<1 

locity of stars, a is the mea n velocity dispersion an d q 
is the axis ratio (|Binnevl[l97i : iBender and Nietdll990l ). If 
(v/a)* > 1.0 then the object is assumed to be rotating. 
We have found that the most massive simulated galaxies 
usually do not rotate or rotate slowly and show signs of tri- 
axiality while less massive galaxies rotate faster and show 
signs of axisymmetry. Thi s statement is in agreement with 
observational studies (e.g. ICappellari" et al.1 (|2007l ) and ref- 
erences therein). Moreover, the majority of rotating galax- 
ies appears to be oblate, rotating around the short axis. 
So for the oblate galaxies observed along the rotation axis 
(and as a consequence seen in a projection with the axis 
ratio q close to unity) the method gives underestimated val- 
ues of the circular speed. It should be noted that when 
observing the rotating galaxies along long axes the circu- 
lar s peed estimate i s slig htly biased towards overestima- 
tion (|Thomas et al.l (|2007 ) reached the similar conlusion). 
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The average deviation for the subsample of oblate galaxies 
seen perpendicular to the rotation axis is biased high by 
Aopt ~ (4.4 ± 1.4)% with RMS = 6.3%. 

To investigate possible projection effects on the results 
of our analysis we have picked one rotating galaxy (the virial 
halo mass is ~ 2.2 x IO^^Mq/i"^) and calculated the surface 
brightness and the velocity dispersion profiles for different 
lines of sight. While the light profiles are quite similar, the 
velocity dispersion profiles may differ significanly when the 
line of sight is parallel to the rotation axis and perpendicular 
to it. We have calculated the average value of the circular 
speed estimates taking into account the probability of ob- 
serving the galaxy at different angles. For the selected galaxy 
the average deviation from the true Vc is about —4.9% and 
the maximum deviation (when observing along the rotation 
axis) is about —25%. 

It should be mentioned that the method under con- 
sideration was designed for recovering the circular speed in 
massive elliptical galaxies and it does not pretend to give 
accurate results for low-mass galaxies. In addition, not so 
many ellipt ical galaxies with a < 150 — 200 km s~^ are ob- 
served (e.g. iBernardi et al]|201Gl ). 

It is convenient to distinguish low and high mass simu- 
lated galaxies by the value of the projected velocity disper- 
sion at the effective radii. Let us call 'massive' galaxies with 
f"(-Rcff) > 150 km s^^. If we apply our analysis to the sub- 
sample of massive galaxies and exclude merging and oblate 
galaxies seen along the rotation axis (the subsample 'MG'), 
we get an unbiased distribution with RMS = 5.4%. The 
resulting histogram is shown in Figure |8] left image, panel 

(A) . Estimations at other radii give slightly more biased and 
slightly less accurate results (Figure [S] left image, panels 

(B) -(D)). 

Thereby we have marked out four subsamples - the 
whole sample without exceptions ('A' - all), the sample with- 
out merging or oblate galaxies seen along the rotation axis 
('G' - good), the subsample of massive galaxies ('M' - mas- 
sive) with CT(_Roff) > 150 km s"'^ and, finally, the subsample 
of massive galaxies when merging and oblate galaxies ob- 
served along the rotation axis are excluded ('MG' - massive 
and good). 

In case of missing or unreliable data on the l ine-of-sight 
velocity dispersion profile [Churazov et all (|2010l ) suggest to 
apply a simplified version of the aforementioned analysis 
(equations ((6}). By neglecting terms 7 and 5 we assume 
that the projected velocity dispersion profile is flat. Then the 
radius at which I{R) oc is the sweet point. The resulting 
histograms for the subsample 'MG' are shown in Figure |8l 
right panel. It can be seen that data on the projected velocity 
dispersion plays noticable role in the analysis if the required 
accuracy is of order of several %. Neglecting its derivatives 
leads to a bias towards underestimated values of Vc {Aopt ~ 
(—4.0 ± 1.1)% at the sweet point) and broader wings/tails 
(RMS = 6.4% at -Rswoct) compared to Figure [S] left panel. 
Nonetheless, if only the surface brightness profile and some 
data on the projected velocity dispersion are available the 
simplified version of the method seems to be a good choice. 



4.2 Simulated galaxies at high redshifts 

We have also tested the same procedure for galaxies at 
higher redshifts. Namely, at z = 1 and z = 2. The frac- 
tion of merging galaxies in the sample is larger at high red- 
shift than at z = and the number of stars in each halo is 
considerably smaller. Nevertheless, results are quite encour- 
aging. For the subsample 'MG' the average deviation of the 
circular speed for the isotropic distribution of orbits at the 
sweet point (estimated via equations Q) from the true one 
is close to zero and the scatter is modest. At redshift z = 1 
the average deviation is Aopt = (—0.3 ± 1.1)% and RMS = 
6.0 %, at z = 2 = (0.9 ± 2.2)% and RMS = 8.0 %. 



4.3 Mass from integrated properties 

Asssuming the logarithmic form of the gravitational poten- 
tial $(r-) — V^lnr + const we can estimate the potential 
$ over some range of radii up to a constant. To calculate 
the potential one needs to know the circular velocity profile. 
If we assume that this profile roughly coincides with the 
isotropic profile Vc^° over some range of radii (let us choose 
R G [0.5-Roff, 3-Roff] as a range of radii easily available for 
observations), we can define: 

<I>°''*= / i-^rfr, (8) 

0.5/?ott 

where R G (0.5-Roff, 3-Roff) and Vc°° can be found from the 
full version of the analysis (the first equation of @ ) or from 
the simplified version (the first equation of As the true 
potential is known we can write cl>*''"° = k . <I>°p* 4- const. In 
the ideal case k = A$*'""'^/A<I>°''* = 1.0. The accuracy of 
such an approach is illustrated in Figure [5] In cyan is shown 
the distribution of subsample 'MG' of galaxies as a func- 
tion of A# = (1 — k) • 100%. Just to remind this subsample 
consists of the massive galaxies with a(J?cff ) > 150 km s~^ 
and merging galaxies as well as oblate galaxies seen along 
their rotation axes are excluded. In case of the full analysis 
the distribution is almost unbiased (the average value of k is 
1.02 ±0.02) with RMS = 10.9%. For the simplified formula 
of we see some offset Tc = 1.09±0.02 and RMS = 11.8%. 
In approximation of small deviations RMS defined for the 
potential calculations is twice as large as RMS for the cir- 
cular velocity calculations because the potential $ scales as 
Vc ■ To compare this approach with previous results let us in- 
troduce the deviation Aopt = i[V^"°f - [K*™]^)/ [V^"""]^ 
estimated at the sweet point -Rswcct. Resulting distribution 
for the same subsample is shown in black in Figure [51 As 
expected the width of this distribution is nearly two times 
larger than the distribution of circular velocity estimates 
(Figure [HJ. 

As we see the gravitational potential can be estimated 
via 14"^° with reasonable accuracy. This fact is in agreement 
with the aforementioned statement that most massive galax- 
ies in the sample have almost fiat circular velocity profiles 
in broad range of radii. 



© 2011 RAS, MNRAS OOO.ITIfTSl 



Simple recipe for estimating galaxy masses 9 



Full analysis, subsample 'UG' Simplified analysis, subsample 'MG' 




deviation A„p^=(V^-Vt™•)/V^'■"^ % deviation A„p^=(V^-Vt™•)/V^'■"^ % 



Figure 8. Left: Distribution of galaxies from the subsample 'MG' (massive galaxies with a-{R^ff) > 150 km s^^ wrhen merging and oblate 
galaxies observed along the rotation axis are excluded) according to their deviations. Deviations are calculated at iJswoct (panel (A)), 
iJcft (panel (B)), 0.5-Rofi (panel (C)) and 2iJ(,ff (panel (D)). Right: The same histograms but for the simplified version of the analysis 
(equations 




Figure 9. Accuracy of the derived potential of massive galaxies (merging and oblate objects seen along the rotation axis are excluded). 
In cyan shown the histogram for the quantity Aj, = (1 — k) • 100%, where k = A^*'"^" / A^°^^ . In black shown the histogram for the 
deviation Aopt of the estimated at the sweet point [V^^^°] from the true one 

|-y^truej2 j^eft: HistO grams for the full version of the analysis 
(equations Q). The average value of k is k = 1.02 ± 0.02 and RMS = 10.3%. The average value of the deviation Aopt is (—0.2 ± 1.9)% 
and RMS = 11.3%. Right: Histograms for the simplified version of the analysis (equations l(6j). k = 1.09 ± 0.02 and RMS = 11.8%, 
Aopt = (-7.2 ± 2.1)% and RMS = 12.7%. 
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Aperture dispersion, subsample 'MG' 
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Figure 10. Distribution of galaxies from the subsample 'MG' (galaxies with (j(R^fi) > 150 km s ^ when merging and oblate galaxies 
seen along the rotation axis are excluded). 



4.4 Circular speed derived from the projected 
dispersion in a fixed aperture 

When data on the velocity dispersion are available only in 
the form of aperture dispersions one can estimate the circu- 
lar speed using the simple relation 



(9) 



where (Jap(< 7?) is the velocity dispersion measured 
within some aperture R. To test this way of evaluating the 
circular speed we have computed the luminosity-weighted 
dispersion within the aperture of radius R (excluding the 
region R < 2Rsoft = 3 kpc where softetning may affect 
results of the analysis) 



0-ap(< R) = 



J I{x)a^{x)xdx 



J I{x)xdx 



(10) 



and calculated the deviation from the true circular 
speed at different radii: at the sweet point -Rswoct for the full 
version of the analysis (equations Q), at J? = 0.5RcS , RcB 
and 2i?cff ■ The resulting histograms for the subsample 'MG' 
are presented in Figure [TO] Comparing with circular speed 
estimations at a single radius (in particular, at the sweet 
point) this method gives a biased result AoptiRes) = (1.0 ± 
1.3)% and noticebly larger RMS (at R^s RMS = 7.8%). 



4.5 Circular speed from X-ray data 

Another way to measure the mass of galaxies comes from the 
analysis of the hot X-ray gas in galaxies. By measuring the 



gas number density n and the temperature T profiles from 
X-ray observations one can find the total mass assuming 
that the gas i s in the hydrostatic equilibrium (HE) (e.g. 
iMathewsl Il978l : iForman. Jones, fc Tuckeij 119851 ) . Assuming 
spherical symmetry, the equation of HE can be written 



IdP 
p dr 



dr 



GM 



(11) 



where P = nkT is the gas pressure (n is the gas number 
density), p = ^nipn is the gas density (nip is the proton 
mass), $ is the gravitational potential, Vc is a circular ve- 
locity and Af is the total mass of the galaxy. In simulations 
the mean atomic weight p is assumed to be equal to 0.58. 
Strictly speaking, assuming the HE one neglects possible 
non-thermal contribution to the pressure, which can be due 
to presence of (i) turbulence in the thermal gas, (ii ) cosmic 
rays and magnetic fields (e.g. IChurazov et al.l l 20081 ). 

To estimate deviations from HE, the so called mass 
bias, we took a subsample of the most massive galaxies with 
M > 6.5 ■ IO'^^Mq. X-ray properties of low mass galaxies 
in the sample are influenced by gravitational softening in 
the central 3-4 kpc and are strongly dominated by cold and 
dense clumps in center. Moreover, we know from observa- 
tions that only the most massive galaxies h ave massive X- 
ray atmospheres (e.g. lO'SuUivan et al.|[200ll ). 

The typical proflles of the gas density and temperature 
extracted from simulations are shown in Figure [TT] We used 
the median value of the electron density Ue and T deter- 
mined in each spherical shell, so that we are free of cold dense 
clumps contamination (jZhuravleva et al]|201ll ). Calculated 
pressure and circular velocity (equation llip are also shown 
in Figure [TT] The spurious feature of simulations is that in 
the central 3-4 kpc cold and dense clumps are strongly dom- 
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Figure 12. Distribution of galaxies according to their deviations 
of estimated circular speed from the true value at i?efF (shown 
in cyan) and at 2RefF (shown in black). Circular speed is derived 
using the hydrostatic equilibrium equation for the hot gas. The 
sample consists of 12 galaxies. 



inating. Even using the median value does not remove these 
clumps, causing strong increase of density and drop of tem- 
perature in the center. These clumps are moving ballistically 
and are not in the HE. 

Deviations from HE Ax = , ,^ " were calculated 

T/truc 
' c 

at RcS (cyan histogram in Figure I12p and 2RcS (black 
histogram in Figure \12\ . The average over the subsample 
value of the deviation at R^s is Ax = (—3.0 ± 1.3)% and 
RMS=4.4%. At 2i?,ff A7 = (-4.0 ± 1.1)%, RMS = 3.8 %. 
V? - 

' True ' X 



The average value of 



V? 

true 



over iicft < R < 2RcH is 



6.8 %. 

To calculate averaged ratio of kinetic energy and ther- 
mal energy — on RcH < R < 2i?efi one should exclude 

-C'therm 

cold dense clumps since their contribution to the kinetic en- 
ergy can be significant. The procedu re to exclude clumps is 
described in lZhuravleva et al.l ()201ll ). In brief, in each radial 
shell, we exclude particles with density exceeding the me- 
dian value by more that 2 standard deviations. An example 
of initial and diffuse projected densities is shown in Figure 



1131 Calculated mean ratio of 



-Bkin 



for diffuse component 



-^therm 

is 4.4 %, which is close to the bias in mass from HE. 



5 DISCUSSION 

In Table [T] we summarize the bias and accuracy of all meth- 
ods discussed above. The sample of simulated galaxies was 
divided into 4 subsamples: the whole sample without ex- 



ceptions ('A'), the subsample 'G'('good') for which merg- 
ing and oblate galaxies observed along the rotation axis 
are excluded, the subsample 'M' of massive galaxies with 
cr(7?eff) > 150 km s~^, and the subsample 'MG' of massive 
galaxies when merging and oblate galaxies seen along the 
rotation axis are excluded. For estimations of the potential 
the bias and the RMS are nearly twice large as those for the 
circular speed estimations. To avoid possible confusion all 
values in the table are associated with Vc - estimations. 

In case of the subsample 'MG' the estimation of the 
circular speed at the sweet point with help of equations ((4]) 
gives the unbiased result (Aopt — 0%) and reasonable accu- 
racy (RMS ~ 5 — 6%). To test whether the unbiased average 
is not just a coincidence we have performed a 'Jack knife' 
test. The resulting average for randomly chosen subsamples 
is less than 1%. The subsample 'MG' consists of 106 objects 
and the statistical uncertainty in this case is about 0.9%. 

For the subsample 'M' of massive galaxies (127 objects, 
26 of them (13.3%) are oblate, 3 of them (2.4%) are with 
ongoing merger activity) we also got almost the unbiased 
average (Aopt ~ (0.2 ± 1.2)%). From an observational point 
of view merging objects can be easily excluded while infor- 
mation on the 'oblateness' of galaxies may not be available. 
If we exclude merging galaxies from the subsample 'M' we 
get the average value of the deviation Aopt ~ —0.7%, RMS 
= 5.9%. So the result is almost unbiased. But if run the 
'Jack knife' tests we get on average slightly underestimated 
values of the circular speed with |Aopt| less than 1.5 %. 

The method is not restricted to nearby galaxies, it also 
allows to recover the circular speed for high-redshift ellipti- 
cals. The circular speed estimate averaged over the subsam- 
ple of massive and slowly or non-rotating simulated objests 
(mergers are excluded) at z = 1 is Aopt = (—0.3 ± 1.1)% 
with RMS — 6.0 %, at 2: = 2 the average deviation is 
A^ = (0.9 ± 2.2)% and RMS = 8.0 %. So the K-estimates 
are also almost unbiased with modest scatter of 6 — 8% as 
in case of subsample 'MG' at z = 0. 

While derivation of equations and ^ is based on 
the assumption of the logarithmic form of the gravitational 
potential we have shown that the circular speed estimate 
at the sweet point is still reasonable even if true circular 
velocity is not flat. 

The case of slowly changing Vc with radius can be il- 
lustrated by the following example. If we assume that Vc 
varies with radius as a power law along with other quanti- 
ties I{R) oc (Tp oc R~^ , P = const we end up with the 
following relation between Vc and ap (from Jeans equation): 



V!{R) 



aliR) 



1 -f g -f 7 - 2/3 r[fi-r[ 



l + Q + 7 1 



(1-/?. 



a + 7 
l+a+7 



) r[i±^^].r[^] 



(12) 



where r[x] is the gamma function. This relation is insensitive 
to the anisotropy parameter /3 when a+'y = 2. One can hope 
therefore that for slopes slowly varying with radius the sweet 
point will be located at the radius where this condition is 
met. Substituting q + 7 = 2 in equation (|12p yields the 
relation between Vc and ap which coincides with equation 
((5| for isotropic orbits for a = 2. Deviations of a from 2 by 
10% cause modest ~ 3% variations in Vc- 

Simulated galaxies are of course more complicated 
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Figure 11. Profiles of hot gas electron number density, temperature, pressure and circular velocity of simulated galaxy. Dotted vertical 
curves show the upper and lower limits on R. Vc plot: solid and dashed curves show mass from HE and total mass from simulations 
respectively. 




Figure 13. Projected number density of hot gas in simulated galaxy. Left: initial density, right: density with excluded clumps. 
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Table 1. Summary of the methods discussed. 
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Figure 14. Distribution of high-redshift galaxies from the subsample 'MG' (massive galaxies with (T(Rcfr) > 150 km when merging 
and oblate galaxies observed along the rotation axis are excluded) according to their circular speed deviations. Deviations are calculated 
at iJaweet (panels (A), (C)) and at i?,cft (panels (B), (D)). 



than the above example. For our sample we have inves- 
tigated possible correlations between the deviation Aopt 
of the estimated Vc from the true one and local (at 
iiswcct) slopes of the velocity dlnT4''"°/dlnr, surface bright- 
ness a = —d\nI{R)/d\nR and velocity dispersion 7 — 
—dlna^/dlnR profiles. There is no obvious correlation be- 
tween Aopt and a or 7. We do see a weak linear trend in 
Aopt and dlnVc'^^'^ /dlnr, although it is much smaller than 
the scatter in Aopt- Most of the galaxies in the sample Vc{R) 
slowly declines with radius near Rsvaet (see Figure [T|. How- 
ever, even after subtracting this trend, the RMS-scatter in 
Aopt is reduced from 5.4% to 5.0%, i.e. only by 0.4%. 

Comparable results are obtained using J [V^^°] ^ jr dr 
over [0.5iicff , 3i?cff] as an estimator of the gravitational po- 
tential. 



The simplified version of the analysis (equations ([H])) 
at the sweet point gives almost the same result as at the 
effective radius. So if one has no enough data to calculate 
all necessary for applying equations Q derivatives it makes 
sense to derive from the first formula of ((Gjl and use 
V^°{Rqs) as an estimation of the circular speed. The qual- 
ity of such approach depends on the 'quality' of the sample. 
In case of non-interacting and almost spherical galaxies the 
RMS is about 7% and the bias is about (-4 ± 1.1)%. As- 
suming flat projected velocity dispersion profile leads to the 
underestimation of the circular speed. If data on the line-of- 
sight velocity dispersion allow to estimate the overall trend 
Aa / AR it may reduce the bias. 

In general we can expect the sweet point to be not far 
from the radius R2 where —d\nI{R)/d\nR ~ 2. Indeed for 
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a smooth surface brightness profile which gradually steepers 
with radius an integral j I{R)RdR diverges at low or high 
limits for —din I {R) /din R greater or lower than 2 respec- 
tively. Therefore one can hope that at the radius R2 the 
contributions to the integral of R < R2 and R > R2 to be 
comparable and R2 ~ ReS- Thus, -Rsweet ^ ReS - E-K- for a 
Sersic model with index n (|Graham and Driverl (120051 )1 

dlnljR) ^_^J R 
dlnR \Rcti j 

and it can be easily seen that the sweet point for the 
circular speed estimation is of order of the effe c tive r adius. 
Moreover, as it was shown in [Ch urazov et al.l lj2010h (Ta- 
ble 4) for Sersic models the stellar anisotropy is close to 
minimal at about 0.5i?cff and this radius can be used as the 
sweet point for the circular speed determination. 

We have tested the statement that -Rswcct ~ -R2 ~ RcS 
on the sample of the simulated objects. If the slope of the 
surface brightness profile is close to —2 over some range of 
radii or a = — din I (R)/ dlnR is not monotonic then there 
is an ambiguity in selecting -Rsweet and R2- To avoid this 
ambiguity we have smoothed I{R) and <t{R) using the width 
of the window function A/ = — 1.0. As a result a{R) 
has become monotonic for majority of objects and newly 
determined -Rswoet , R2 follow the relationship -Rswcct ~ -R2 ~ 
Reff. However, a significant smoothing of data leads to a bias 
in estimating the circular speed Aopt ^ —2% at both Tiaweet 
and i?2. 



6 CONCLUSIONS 

Being an important issue, the total mass estimation for ellip- 
tical galaxies is often quite difficult, especially for galaxies 
at high redshift. We used a large sample of cosmological 
zoom simulations of individual galaxies to test a simple and 
robust procedure (see equations Q, (|6}) based on the sur- 
face brightness and velocity dispersion profiles to estimate 
the circular speed and therefore the total mass of a mas- 
sive galaxy. The method is very simple and it does not re- 
quire any assumptions on the stellar anisotropy profile. For 
massive ellipticals without significant rotation at redshifts 
z = — 2 it gives an unbiased estimate of the circular speed 
(the bias Aopt(iiswect) is less than 1%) with 5-6% scatter. 
Therefore this method is suitable for the analyze of large 
samples of galaxies with limited observational data at low 
and high redshifts. The method works best for the most mas- 
sive ellipticals {a{Rcff) > 200 km s~^), which in the present 
simulations have almost isothermal circular velocity profiles 
over broad range of radii. 

The method should be applied with caution to merging 
galaxies where the circular speed can be significantly over- 
estimated. For rotating galaxies seen along the rotation axis 
the procedure gives substantially underestimated 14. 

The best estimate of the circular speed is obtained at 
a sweet point 7?aweet where the sensitivity of the recovered 
circular speed to the stellar anisotropy is expected to be 
minimal (see section [2]). The -Rsweet is expected to be not 
far from the projected radius where the surface brightness 



declines approximately as / oc R~^. This radius is in turn 
close (within factor of 2) to the effective radius Rett of the 
galaxy. Our tests have shown that the accuracy (RMS scat- 
ter) of the circular speed estimates at 0.5 — 2 Rett is 5 — 7% 
for most massive ellipticals. 

An even simpler method - based on the aperture ve- 
locity dispersion (equations ((Oj, H10[l) - is found to be less 
accurate, although the results are still reasonable. For ex- 
ample, for massive galaxies without significant rotation the 
sample averaged deviation of the circular speed at the effec- 
tive radius is AoptiRcfi) = (1.0 ± 1.3)% with RMS ~ 8%. 
Other fiavors of the circular speed estimates are described 
in Section [Til 

Using the same simulated set we have also tested the 
accuracy of the circular speed estimate from the hydrostatic 
equilibrium equation for the hot gas in massive ellipticals. 
We found a negative bias at the level of 3 — 4% and the scat- 
ter of ~ 5%. The presence of bias is caused by the residual 
gas motions. 

Given the simplicity of the described method (Section 
the low bias and modest scatter in the recovered value 
of the circular speed, it is suitable for the analysis of large 
samples of massive elliptical galaxies at low as well as high 
reshifts. 
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